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Abstract 

It is shown that layered superconductors are subjected to a phase transition 

at zero temperature provided the order parameter (OP) reverses its sign on 

the Fermi-surface but its angular average is finite. The transition is regulated 

by an elastic impurity scattering rate 1/r. The excitation energy spectrum, 

being gapless at the low level of scattering, develops a gap as soon as the 

scattering rate exceeds some critical value of 1/t*. 
PACS numbers: 74.20.Fg, 74.25.Jb, 74.62.-c, 74.72.Bk 
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It has been shown recently that non-magnetic impurities may destroy the d-w&ve super- 
conductivity in the same manner as magnetic impurities do it in the s-wave conventional 
superconductors |I|-£§. However, this pure <i-wave scenario is not supported by the existent 
experimental data, at least in YBaCuO samples, subjected to radiation damage or doping 
0. One might anticipate that T c should drop abruptly above a certain critical concentration 
of the defects. Instead, critical temperature has been found just to decrease gradually in the 
whole range of a residual resistivity variation. The above discrepancy has been considered 
as a compelling counter-argument against the whole idea of exotic pairing in this compound. 
However, in this note we make the case that the above mentioned experimental results do 
not in fact contradict the alleged non-trivial structure of the order parameter (OP), provided 
the latter is determined by the symmetry of the crystal and, hence, deviates from the exact 
<i-wave form. In this situation, as will be demonstrated below, an initially gapless excitation 
spectrum of a clean superconductor may acquire a gap due to a scattering of electrons by 
non- magnetic impurities. An elastic scattering gives rise to this paradoxical effect, once 
the OP in the clean material possesses the nodes on the Fermi surface whereas its angular 
average still does not vanish, and the scattering rate exceeds some critical value. 

The results of recent Josephson tunneling experiments [Q and the measurement of critical 
temperature versus residual resistivity variation || may be interpreted as an evidence in 
favor of this very structure of the OP. Indeed, the Josephson tunneling data in the corner 
geometry imply the sign reversal of the OP on the Fermi surface. On the other hand, as is 
already mentioned, the growth of the residual resistivity is accompanied by a slow decrease 
of critical temperature. This type of behavior rules out the nullification of the OP averaged 
over the Fermi surface as we shall show later. 

It should be noted that recent measurements of a Josephson current in c-direction by Sun 
et al. do not conform easily with the phase alternating OP. The values of the net current 
found in this experiment for the heavily twinned samples in plain geometry are several orders 
of magnitude higher than theoretically estimated. The interpretation of these results may 
prove to be a subtle matter, however. A possible source of the finite Josephson current could 
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be a symmetry violation caused by surface defects and the surface itself together with the 
phase self-adjustment. For other plausible scenarios this effect see preprints J7|. For now we 
assume that the OP sign reversal, as well as the non-zero value of (A), is reliably established 
in the experiments cited above @||. 

We analyze the problem in the framework of the anisotropic BCS model. The OP is 
determined by self-consistency equation: 

A W ^E/^fl/- W f (1) 



Here the summation goes over the Matsubara imaginary frequencies whereas the integration 
is restricted to the Fermi surface, which is presumed to be a cylinder. The kernel V(<f>, <p') 
describes the interaction of electrons on the Fermi surface. The renormalized Matsubara 
frequency fj n and the frequency- dependent OP A n (0) are related to the corresponding bare 
quantities r\ n = (2 n + 1) 7r T and A(0) via Abrikosov-Gor'kov (AG) equations: 

Vn-^( , = )=Vn (2) 



*»(*)--( , Ara(0/) ) = m (3) 

The angular brackets in the AG equations denote the angular average. For the sake of 
simplicity we consider an isotropic scattering only. For the same reason we have disregarded 
any dependence of the OP and the electronic interaction potential on the c-component of 
momentum. 

Let us first examine the behavior of the critical temperature with respect to variation of 
the scattering rate. In the vicinity of critical line the equations (|TJ-^J) may be linearized and 
the infinitesimal order parameter may be eliminated. The result reads: 

1 = g(T,r) ((l - fiT^Vy'V^)) (4) 

Here we defined a linear operator V related to the electronic interaction kernel: 

vm = ( (5) 

J 271 
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This operator maps 1 onto the function V(4>) 



V(<P) = Vfatfl-f- (6) 

J 271 

The functions f(T, r) and g(T, r) in eqn. (|J) are defined by the following expressions: 

The summation in the expression for f(T, r) is limited by an ultraviolet cutoff: r] n < e. The 
latter is related to the critical temperature T c q of the clean superconductor as follows: 

r c0 = ^exp(-7r/Fo) (9) 

7T 

where Vq is the maximal eigenvalue of the operator V and 7 is the Euler constant. In the 
dirty limit (tT c0 <C 1) the approximate solution of eqn. (f|) reads 

T c = 7 e( 7 er/2)^ 1 exp(-7r/(y(0))) (10) 

where k = (V 2 )/(V) 2 . Within the same approximation A(0) oc V(^). Hence, the exponent 
k may be expressed in terms of the OP: k = (A 2 )/(A) 2 . The power law (|TD|) has been 
derived in the sixties by P. Hohenberg || under the assumption of a weak anisotropy. We 
have found that it is the high scattering level that really matters. Notice that the exponent 
k becomes infinite when (A) — > 0. A more scrupulous analysis [0-0 shows that if (A) = 
the critical temperature vanishes at some finite scattering rate proportional to the critical 
temperature of the clean superconductor: 2tit c T c0 = 1. 

Equation ( |T0D might be interpreted directly as the relation between the critical tempera- 
ture and the residual resistivity: T c (p) oc p K_1 , provided that the effective number of carriers 
does not depend substantially on the concentration of defects. 

Let us study equations at temperature equal to zero. In this case the summation 
in eqn. ([I]) should be replaced by integration over the continuous variable 77 and the domain 
of functions fj(r]) and A(0, rj) extends into the whole positive half-axis of the same variable. 



We are interested in the behavior of the above functions at 77 — > 0. The values 7/(0) and 
A(0, 0) determine the density of states (DOS) on the Fermi surface, vanishing if 77(0) = 0. 

We are going to prove the following three statements concerning the DOS on the Fermi 
surface: 

1. Let (A(0)) 7^ and A(0) has no nodes. Then 77(0) = for any r. 

2. Let (A((f))} = 0. Then fj(0) > for any r > (271^) -1 . 

3. Let (A(0)) 7^ but A(0) possesses nodes. Then 77(0) = for r < t* but 7/(0) > for 
r > r*. The equation for T* as a functional of A(0) will be derived below. 

Before proceeding to the proof let us remark that according to eqn. (|2|) the following 
separation of variables takes place: 

A(0,77)=A(0)+a(77) (11) 

and the AG equations may be rewritten in terms of the function 17(77) j us ^ defined: 

77(77) (l - - / 1 \ I = 77 (12) 

V T \y/(A(<l>)+<T( V ))* + ri( V )*/j M^/(A(0) + a(77))2 + 7/(77)2/ 
The first proposition stems from eqns. (|T2]JT3|) straightforwardly. Really, let 77(0) 7^ 0, then 



the expression in the brackets vanishes and, consequently, the r.h.s of eqn. (P^ ) must have 
a root at 77 = 0. The latter is impossible, however, if A(0) does not reverse its sign. Hence, 
if A(0) has no nodes 7/(0) = 0. 

Let us skip the second statement for a moment and consider the case when the OP in a 
clean superconductor does reverse its sign but its angular average is finite. 

We study, first, the limit of small impurity concentration (r(A(0)) ^> 1). Collecting 
terms proportional to 1/r in eqn. (]T3|) one can find out that a = 0(1/t). In other words, 
renormalization of the order parameter is small. In particular, its nodes remain almost at 
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their original locations. However, renormalization of 77(77) is not small due to the logarithmic 
divergence of the term, proportional to 1/r in the l.h.s. of eqn. (|T^) at 77(77) = 0. A closer 
examination of the l.h.s of eqn. (|T2|), considered as the function 77(77) (inverse to 77(77)), shows 
that it departs from the coordinate origin with an infinite negative derivative and, after 
reaching its minimum, crosses the abscissa axis at 770 = A' L exp (— rAj/2), where 1/Aj = 
(1/ I A; I +1/ I A' 2 |), whereas log A' L = (| A[ | log | A[ | + | A> 2 | log | A' 2 |)/(| A[ | + | A 2 |) 
and A[ 2 denote the derivatives of the order parameter at its nodes. Among the two roots, 
only 770 resides on the physical sheet |J. The first part of the statement 3 is proved. 

In the opposite limit of large scattering rate (1/r 3> T c0 ) the asymptotic solution of eqns. 
(00) can be found in the following form: 

V(V) = V[1+ 1 1 9 (14) 
V rV<A(0)> 2 +W 

^ = . m (15) 

ry/(A(<l>)) 2 + rP 

Plugging this solution into the self-consistency relation at T = we derive the following 
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vrA(0) = log (re) VA(<P) - log (r(A» V(<P)(A) (16) 

Assuming that log re ^> , where Vq is the the maximal eigenvalue of the operator V, as 
previously, its solution can be found explicitly: 

A(0) = (A)V(<P)/(V); (A) = 2r K - 1 e- K exp(- 7 r/(V)) (17) 



From eqns. (|14| , |T5"|) it follows that 77(0) = and a(0) = 1/r 3> A. Since the renormalized 
OP A(0) does not have any nodes, no divergence can arise in eqn. (0). In this way a 
self-consistency of the obtained solution is guaranteed. 

Summarizing, 7/(0) 7^ in a clean superconductor and it vanishes when a scattering rate 
substantially exceeds T c0 . Hence, there should exist some special value of r = r* at which 
77(0) first turns into zero. Since a(0)\ T=T = 1/r*, the value of r* can be found as a functional 
of A(0) by means of the following equation: 



(A(0)/(A(0) + l/r*)> = (18) 

The above equation should be supplemented by the self-consistency and the AG equations 
(|I|,|T^, [13]) with t = r*. Using this ansatz one can find r*, at least in principle, as a functional 
of the OP in a clean superconductor. 

To get more visible results we consider a limiting case (V) 2 <C (V 2 ). If this strong 
inequality is satisfied, the average OP should be small compared to the amplitude of its 
variation and one can expect that the dependence T c (r) is close to that for (A) = 0. It 
means that T c (r) reaches almost zero value as r — > r c , but then has a power-like tail: 
T c {t) ~ t k_1 with k ^> 1. Since the point r* separates the domain of r with the d-like 
behavior from that of the s-like one, it is natural to conjecture that is close to r c . The 



direct calculation justifies this guess. We assume that 1/r* J (A 2 ) ^> (A), then eqn. (|I8D 
yields 

r* = (A)/(A 2 ) (19) 

Plugging the solution of eqn. (]16[) into the expression for (A 2 ) of the above relation one 

gets: 



T*(A> (l0g(T*<A»/7T) 



l--log(r ir e)v) V(0) 

7T / 



1 (20) 



Since t*(A) is assumed to be small, log(r^e)/7r must be close to the reciprocal maximal 
eigenvalue V^ -1 of the operator V in order to satisfy the above equation. Therefore, 
coincides with r c for (A) 2 / (A 2 ) — > 0. This completes the proof of the statements 2 and 3. 

The properties of the AG equations, which have been just established, give us a clue for 
an investigation of the behavior of the angular-dependent DOS. The latter may be found 
via the solutions of the AG equations, analytically continued from the Matsubara sequence, 
defined on the upper half of the imaginary axis of the frequency complex plane, to the real 
axis according to the formula: 

u s (e, 0) = u n (e, ^Ke \ 6 M } (21) 

e 2 (e) - A 2 (6,e' 



Here v s (<f>) and v n {4>) denote the DOS in the superconducting and the normal states respec- 
tively, and the functions e(e) and A(0, e) represent the result of the analytical continuation 
mentioned above. Namely, e(ir] n ) = ifj n , A(4>, irj n ) = A ra (0). At the Fermi-level e(0) = ifj(0) 
for r > r* and e(0) = for r < r*. Hence, zA;(</>, 0) is finite for r > and vanishes identically 
for r < r*. It also vanishes at r = oo. Therefore, it has a maximum at some r > r*. 

It is natural to suppose, that a gap in the excitation spectrum does persist for all r < r*. 
We can confirm this conjecture at least for r <C r*. The solution (|14| , |15| ) of the AG equations 



can be employed in this case with the same substitution rj = —ie and fj = —ie. Plugging 
this solution into eqn. ( f2~l|) we find with the precision up to rA: 

v s (e, 0) = z/ n (e, 0)^e J . i (22) 



2 - (A) 2 



This is exactly the same formula as for an isotropic superconductor with the average value 
(A) playing the role of the isotropic gap. This result agrees with the Anderson isotropisation 
theorem |IT| |. Note, however, that neither h> s (e,<j)) nor the OP become isotropic even in the 



extra dirty limit. 

Evaluation with a higher precision shows, that the standard singularity in DOS 
turns into a finite maximum of the height ~ ((A)r) _1//2 , smeared over the interval 5t ~ rA 2 
and varying with angle. Nevertheless, the threshold character of the DOS dependence on 
energy with the threshold approximately equal to (A) remains unaffected. 

We have shown that the gap vanishes at r = and it also vanishes at r = 0. Therefore, 
the gap reaches its maximum value at some r < r*. 

Thus, we propose to look for a rather peculiar phase transition, which can exist at zero 
temperature and is regulated by impurity concentration. One can try to find it experi- 
mentally, scrutinizing the dependence of quasiparticle tunneling rate or the temperature 
correction to the penetration depth on the residual resistivity in single crystals of YBaCuO, 
doped by Pr or subjected to the radiation damage. A transmutation of the gapless tunneling 
current- voltage characteristics at small resistivity into the curves with a threshold feature at 
higher values of residual resistivity and an analogous conversion of the temperature depen- 
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dence of the penetration depth from the gapless to the activated one would clearly indicate 
the presence of the transition. An ideal experimental setting would include the controlling 
Josephson tunneling measurements in the corner geometry on the same samples to verify 
whether the OP changes its sign. 

In conclusion, we have shown that a phase transition should take place at zero temper- 
ature and a special value of the impurity scattering rate 1/r* in a layered superconductor, 
provided the OP in the clean limit does change its sign on the Fermi-surface but its angular 
average is finite. We have argued that, according to the recent experimental observations, 
the OP with the required properties is likely to be an inherent feature, at least for YBaCuO. 
This transition is characterized by a gap generation in the excitation spectrum for r < r*. 
The energy gap grows from zero at r = to its maximum at some r > and vanishes 
again in the extra dirty limit. The DOS at the Fermi- level turns into zero in the clean case 
and at r = r*, reaching its maximum at an intermediate value of r. 

One of the authors (V.L.P) thanks ETH, Zurich, J. Blatter and especially D. Pescia for 
the hospitality and support extended to him during the time when this work was in progress. 
Thanks are due to H. Monien and D. Geshkenbein for discussion and indication of several 
references. 
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